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ABSTRACT 

We analyze both the early and late time radio and X-ray data of the tidal 
disruption event Swift J1644+57. The data at early times (< 5 days) necessi- 
tates separation of the radio and X-ray emission regions, either spatially or in 
velocity space. This leads us to suggest a two component jet model, in which the 
inner jet is initially relativistic with Lorentz factor F ^ 15, while the outer jet is 
trans- relativistic, with F < 1.2. This model enables a self-consistent interpreta- 
tion of the late time radio data, both in terms of peak frequency and flux. We 
solve the dynamics, radiative cooling and expected radiation from both jet com- 
ponents. We show that while during the first month synchrotron emission from 
the outer jet dominates the radio emission, at later times radiation from ambient 
gas collected by the inner jet dominates. This provides a natural explanation to 
the observed re-brightening, without the need for late time inner engine activity. 
After 100 days, the radio emission peak is in the optically thick regime, leading 
to a decay of both the flux and peak frequency at later times. Our model's pre- 
dictions for the evolution of radio emission in jetted tidal disruption events can 
be tested by future observations. 

Subject headings: black hole physics — galaxies: jets — galaxies: nuclei 

radiation mechanisms:non-thermal 



1. Introduction 

A stray star, when passing near a massive black hole (MBH) can be torn apart by 
gravitational forces, leading to a tidal disruption event (TDE). Such an event would be 
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observed as bright emission from a previously dormant MBH, as it is being f ed by tem- 



porary mass accretion established a fter the tidal disruption of a passing star (IHillsl Il975 



Reeslll988l : lEvans &: Kochaneklll989l ). On 25 March 2011, an unusual transient source Swift 
J1644 49.3+573451 (hereafter Sw J1644+57) was reported, potentially representing such an 
event (IBurrows et al.ll2011t iLevan et al.ll201ll ). This event was found to be in positional co- 
i ncidence (5, 0-2 kpc) with a previously dormant host gala xy nucleus, at redshift z = 0.354 
Levan et aPboili livuchter et al.lboili ISerger et al.lboill ). 



The rapid variability seen in the X-rays, of ~ 78 s ( iBurrows et al.l 1201 ll : iBloom et al. 



201 ll ). implies a compact source size < 0.15 AU, which is a few times the Schwarzschild 



radius of IO^Mq MBH (see also iMiller fc Giiltekinl I20111). When cornbined with the ver y 



high 7-ray and X-ray luminosity, ~ 10 erg s (IBurrows et al.l l2011t iBloom et al.l l201ll ) 



which is 2 - 3 orders of magnitude above the Eddington limit of such a MBH, it was 
concluded th at the X-ray emiss ion must origina te from a relat i vistic jet of Lorentz fac 



Quataert fc Kasen 



2012 



tor r > 10 (iBloom et al.ll201lr) (see, however, iKrohk fc PiranI I2OIII : lOuved et al.l I2OII 



Socratesll2012l . for alternative models). 



While early works that investigated the expected signal (optical/UV emission) from 
such an event were focused on the signal from the accreting material of the stellar de- 
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bris (lReeslll988l: iLoeb fc U 



Strubbe &: Quataert 



jet was considered 



Metzger et al.l 12012 



mei 



19971 : IUlmerlll999l : iBoedanovic et al.l 120041 : iGuillochon et al. 



2009), recent 



Giannios fc Metzgei 



De CoUe et al 



y, the observational signature from a newly formed 



2011 



2012 



van Velzen et al 



2OIII : IWang fc Cheng|l2012 



Stone &: Loebll2012l ). The basic mechanism sug- 



gested in these works is similar to the mechanism that is thought to oper ate in gamma-ray 
bursts (GRBs), namely ener gy dissipation by either internal shock waves ( iRees fc Meszaros 
I994J : iPaczynski fc Xul Il994j ) or by an external (forward) shock wave, accompanied at early 



stages by a reverse shock wave (IRees fc Meszarod Il992l : iMeszaros fc ReesI Il993l : iPiran et al. 



I993I : ISari fc PiranI I1995I ). These shock waves, in turn, are believed to accelerate particles 
and generate strong magnetic fields, thereby producing synchrotron radiation, accompanied 
by synchrotron self Compton (SSC) emission at high energies. 

Indeed, shortly after its discovery, an extensive radio campaign show ed that the X-ray 
emission is accompanied by bright radio emission (IZauderer et al.l 1201 ll ). which was inter- 
preted as synchrotron emission from the jetted material, thereby supporting this hypothesis. 
However, a careful analysis reve aled that the radio em itting material propagates at a more 
modest Lorentz factor, F > 1 ( IZauderer et al.ll201ll ). and therefore the X-ray and radio 
emission cannot have similar origin. This had led to the suggestion that the X-rays may 
originate from internal dissipation ( IWang &: Chengl |2012| ). while the radio may originate 



from the forward shock that propagates into the surrounding material (iMetzger et al 



2012 
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Cao fc Wanell2012h 



Late time radio monitoring, extending up to ~ 216 days (iBerger et al.ll2012l ). revealed 
an unexpected behavior: after about ~ 30 days (observed time), the radio emission showed 
re-brightening, which lasted up to ~ 100 days, after which the radio flux decayed. This 
re-brightening was not accompanied by re-brightening in the X-ray fl ux, a nd is not expected 
in the context of the forward shock models. This led iBerger et al.l (120121 ) to conclude that 
the radio re-brightening resulted from late time energy injection. Thus, a comprehensive 
model that considers the temporal as well as combined radio and X-ray spectral data is still 
lacking. 

Any such model must take into account the fact that the emitting material is mildly 
relativistic at most. First, the radio emitting material propagates at Lorentz factor F > 1. 
Second, while the X-ray emitting material propagates at initial Lorentz factor F ^ 10 — 
20, its velocity becomes trans-relativistic on the relevant time scale of tens of days, as 
surrounding material i s collected. The dynamics of such trans-relativistic propagation was 
recently considered by iPe'erl (120121 ). 



Here, we propose a new model that considers simultaneously the emission of both the 
radio and the X-rays, their spectrum as well as temporal evolution. We re-derive the con- 
straints set by both radio and X-ray observations, and confirm that indeed at early times 
(first few days) these must have a separate origin. We calculate the dynamics of the X- 
ray emitting plasma as it collects material from the surrounding and decelerates. We show 
that after ~ 30 days (observed time), synchrotron emission from this plasma peaks at radio 
frequencies, thereby providing a natural explanation to the re-brightening seen at radio fre- 
quencies at these times, without the need for late time internal engine activity. Moreover, we 
show that the decay of the radio flux after ~ 100 days is naturally explained by synchrotron 
self absorption. At this stage, the flow is in the trans-relativistic regime (F — 1 ~ 1). 

This paper is organized as follows. In ^ we carefully revise the observed data of 
both the radio and X-ray emission at early times (up to five days). While our treatment 
is more general than previous works, we confirm earlier conclusions that indeed the radio 
and X-ray emission must have separate origins. In ^ we consider the temporal evolution 
of the X-ray emitting material as it slows and cools, and show that it can be the source of 
the re-brightening at radio frequencies seen after ~ 30 days. We further consider radiative 
cooling in §3.21 We show that the very late (> 100 days) decay of the radio flux is naturally 
attributed to emission in the optically thick regime: as the electrons cool, eventually the 
peak of the synchrotron becomes obscured. We compare our model to the late time radio 
data of Sw J1644+57 in §U before summarizing our main results in ^ 
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2. Early time radio and X-ray emission and its interpretation 



Sw J1644+57 is a long-lived (duration > months) X-ray outburst source, accompanied 
by bright radio emission, interpreted as synchrotron radiation. During the first few days of 
observations, the iso tropic X-ray luminos ity ranged from ~ 3 x lO^^ergs"^ to a peak as high 
as ~ SxlO'^^ergs"^ ( Burrows et al.|[2011 ) with average X-ray luminosity ~ /ewx 10^^ erg s^^. 
The X-ray emission peaked at frequency ^ 2 x 10^^ Hz, with e stimated uncertainty up to 
two orders of magnitude above this value (jZauderer et al.ll201ll ). During the first few days, 
the X-ra y light curve was cornplex and highly va riable, with variability time scale as short as 
~ 100 s jBurrows et al.lboill : Ihevan et al.lboilh . 



This source triggered a radio campaign, that began a few days after its initial discovery. 
Radio observations showed that the peak frequency occurs at ~ 8x lO^^r] Hz with uncertainty 
1 < < 10, and peak luminosity z/Lj, > lO^^'^ergs"^ during the first few days. The spectral 
energy distribution (SED) at the radio band (< 345 GHz) at At°^ ^ 5 days is well described 
by a power law, F^, oc u^''^ up to -Fy(|iy=345GHz) ~ 35 mJy The steep power law index 
requires self-absorbed synchrotron emis sion, with self absorption frequency z/^ > 10^^ Hz 
(jZauderer et al.ll201ll : iBerger et al.ll2012l ). Within two weeks, the X-rays maintained a more 
steady level, albeit with episodic brightening and fading spanning more than an order of 



magnitude in flux, while th e low frequency (radio) emission decreased markedly (ILevan et al. 
20 111 : IZauderer et al.l[201l[ ). 



2.1. Interpretation of early radio emission 



The radio emission observed from Swift J1644+57 is assumed to have synchrotron origin. 
The synchrotron emitting plasma can be described by five free parameters: the source size 
R, bulk Lorentz factor F, total number of radiating electrons N^, magnetic field strength 
B, and the characteristic electron Lorentz factor (as is m easured in the plasma frame), 7e. 



Calculations of the values of these parani e ters a ppear in IZauderer et al 



(1201 ih . Here, we 



generalize the treatment in IZauderer et al.l (120111 ) by removing the equipartition assumption 
used in that work. 

Existing data provides the following four constraints. The observed ch aracteristic fre- 
quen cy and total luminosity of synchrotron emission, i/^ and uL^ are given by (IRybicki fc Lightman 
1979h 



,ob 



I'm 



3 2]- ^ 4 20 X 10** 5 7,^r = 8 X 10^^ Hz: 



1^ 



f B-ilV = Sx 10" 
syn - iNearc-ilUBV^ = 1.06 x 10-15 N^-fj B^T^ = 3x 10^^ erg s-\ 
Here and below, Pgyn is the total synchrotron radiation power, ay is the Thomson cross 
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section, Ub = -B^/Stt is the magnetic energy density and CGS units are used. In deriving 
equation ([1]), spherical explosion was assumed. 

A third condition is the synchrotron self absorption frequency, Ua ^ 2 x 10^^ Hz. For 
, < Um., as is the case here, the self absorption coefficient scales as oc z/~^/^ (see 



Rybicki fc Lightman 1979 . equation (6.50)), thus Ty/ry^ = a^/ay^ = {y /umY^^^ ■ Here, 
Ty OC ay is the optical depth, Ty^ = Ty{v = Vrn) = cty^R' , and Ft = R/T is the co- 
m oving size of the synchrotro n emitting region. In calculating ay^ we use equation (6.53) 



m 



Rybicki &: LightmanI ( 119791 ) . which assumes power law distribution of electrons above 7e. 
We use power law index p = 2, although the result is found not to be sensitive to neither 
the power law distribution assumption or the power law index; a similar result is obtained if 
the electrons assume a thermal distribution. Since, by definition, t^^ = Ty{v = Va) = the 
synchrotron self absorption frequency Va is 

^ = 2.46 X IQ-^Nl/^R-^'^B-^'^^"^ ~ - . (2) 

Here, the observed value of the self absorption frequency v'^ ^ 2 x 10^^ Hz is used. 

As a fourth conditio n we use the assumption of ballistic expansion of the source, similar 



to 



Zauderer et al.l (120111 ). We consider constant bulk expansion velocity of the source j3 = 
(1 — r~^)^^^ and head-on emission. At redshift z = 0.354, an observed time Af^^ = 5 days 
corresponds to source frame time At°^/ {1 + z) = 3.7 days. Due to relativistic time compres- 
sion, the actual time that the source would have expanded is [At"^/{1 + z)]/{l — /3). Thus, 
the emission radius is related to the observed time by 

At-" Pc 

The four constraints derived in equations ([I]) - ([3]) are insufficient to fully determine 
the values of the five free model parameters. We therefore choose the source size R at 
At°^ = 5 days as a free variable, and determine the values of the other four free parameters. 
The results of our calculation are shown in Figure [TJ We present the bulk Lorentz factor F 
(black), bulk momentum F/3 (cyan), magnetic field B (green), characteristic electron Lorentz 
factor 7e (red) and the number density of the radiating particles, Ue (blue). 

In order to constrain the allowed parameter space region, we use two additional assump- 
tions: (1) In Fermi-type acceleration, the typical Lorentz factor of the energetic electrons 
7e < ^p/fTie = 1836 in the rest frame of the plasma; (2) As the emission radius is large, 
the magnetic field must be produced at the shock front. Thus, the ratio of magnetic energy 
density, B"^ /9>-k to the photon energy density, L/(47ri?^F^c) (often denoted by e^), is smaller 
than unity. These constraints are shown by the dashed lines in Figure [TJ 



-2 I ^ I , ■ 

10^' 10' 



R [cm] 



Fig. 1. — Dependence of the free model parameters: the bulk Lorentz factor, F (black), 
the bulk momentum, F/3 (cyan), the magnetic field, B (green), the typical Lorentz factor of 
electrons, 7e (red) and the number density of electrons, rie (blue) on the emission radius of the 
radio photons, R, under the assumption that the radio photons originating from synchrotron 
emission. 



After adding these two constraints, we conclude that the radio emission zone at early 
times fulfills the following conditions: (i) The emission radius is in the range 3.0 x 10^^ cm < 
i? < 1.0 X 10^^ cm. (ii) The outflow is trans-relativistic, with 1.04 < F < 1.2, and 0.24 < 
f3 < 0.55. (Hi) The magnetic field is poorly constrained by the data, and could range from 
5.0 X 10^^ G to 7.0 G. (iv) The electrons are hot, with minimum random Lorentz factor 
exceeding 7e > 150. (v) The emitting region is dense: 3.0 x 10^ cm~'^ < < 3.0 x 10^ cm~^. 
This result likely excludes shocked external material as the source of the radio emission, as 
the typical densities (even after compressed by mild-relativistic shock waves) are much lower 
than these values, (vi) The total number of radiating particles is at the range 1.0 x 10^^ < 
< 3.4 X 10^^, with likely value of ~ 1 0^^. We find that our results are consistent with the 
results derived by IZauderer et al.l (1201 ll ). 
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2.2. Interpretation of early X-ray emission 



Our underlying assumption is tliat tlie origin of tlie X-ray emission is inverse Compton 
(IC) scattering of the synchrotron radio photons by relativistic electrons. As we show here, 
these electrons must be l ocated in a di f ferent region than the radio-emitting electrons. Fol- 



lowing the treatment by iPe'er &: Loebl (120121 ). two constraints can be put by the data: (I) 



the ratio of the IC and synchrotron peak frequencies is given by 



3 ' ' ' ' 8 X 10" 



and, (II) the ratio of IC to synchrotron peak fluxes is 

-^!/,IC _ {l-'Fy^peak,Ic/^peak,lc) 



n, 



e[ic]rcrT ^ 4 X 10"^ (5) 



Here, 7e[/c] and ne[ic] are the typical Lorentz factor and number density of electrons that 
emit the IC photons, and r and F are the typical size and the bulk Lorentz factor of the 
IC emission region. In estimating the ratio of the peak frequencies and fluxes, we used 

<ak,IC - 2 X 10^8 Hz, i/ob^j^^^y^ ^ 8 X 10" Hz, iyF,^peak,IC ~ 3 X 10^^ erg S'^, and iyF,^peak,syn ~ 

3 X 10^3 g-i^ 

The rapid variability observed in X-rays on a time scale St°^ ~ 100 s constraints the 
size of the emitting region, r. While the time during which this variability is observed does 
not correspond directly to flve days, we use it here as an order of magnitude estimate. This 
variability implies a relation between the emission radius of the X-rays and the bulk Lorentz 
factor, 

^ ' - 2T^c- (6) 



/3 l + z 



Equations (jlj), ([5]) and ([6]) exhibit three restrictive conditions provided by the observed 
data. As there are four free model parameters, F, 7e[/q, ne[ic] (or Ne) and r, a full solution 
cannot be obtained. However, we can apply a similar method to the one used in §2.H namely 
take r as a free parameter and obtain the values of the other three unknowns. The results 
of this analysis are presented in Figure [2l In this flgure, we present the values of F (black), 
le[ic] (red) and ne[ic] (blue) as a function of r, where 5t°'' = 100 s is considered. 

From the results of Figure [2] one can put several constraints on the X-ray emission zone: 
First, there is no strong constraint on the size of the emitting region, r. Any value in the 
range 10^^ cm < r < 10^^ cm, which is compatible with the size of the emission zone of the 
synchrotron photons, is acceptable. The main constraint originate from equation (l6l), as large 
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r [cm] 



Fig. 2. — Dependence of the free model parameters of the X-ray emitting region: the bulk 
Lorentz factor, F (black) the typical Lorentz factor of electrons, 7e (red) and the number 
density of electrons, tie (blue) on the radius of source, r. We assume that IC scattering of 
the radio photons are the main source of X-rays. 

radius implies large bulk Lorentz factor, as F ~ i?^/^/(2.1 x 10^). Thus, value of r ^ 10^^ 
imply F ^ 50, much larger than the value obtained for the radio emission zone. Even value of 
r = 10^^ cm implies F ^ 15, inconsistent with the finding for the radio emitting zone. Thus, 
while the size of the radio and X-ray emitting regions can be comparable, the X-ray emitting 
region must propagate at much larger Lorentz factor than the radio emitting region. Second, 
the IC emitting electrons are not as hot as the radio emitting electrons. For r > 10^^ cm, 
F > 15, which, using equation (jlj) imply 7e[/c] ^ 100. Third, the number density of the 
radiating electrons, 10*^ cm~^ ^ ^ 10^ cm~^ is comparable to the number density of the 
radio emitting particles, and is thus likely too high to be explained by compression of the 
external material. A similar conclusion holds for the total number of the radiating particles. 

Thus, we conclude that while the two emission regions can have a comparable size, the 
X-ray emission zone must have a much larger Lorentz factor than the radio emission zone; 
moreover, the electrons in the X-ray emitting region must be colder than the electrons that 
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emit at radio frequencies. Alternatively, the bulk motion may be similar, but only if the 
X-ray emitting region is at much smaller distance, lO^^ ^ qjxi. These results are summarized 
in Table [T] 



Table 1: Summary of key properties of both radio and X-ray emission zones of Swift 
J164449.3+573451 





R [cm] 


r 




7e 


Tif. [cm ^] 


K 


Radio emission 














zone 


3 X 10^5 - 1016 


<1.2 


150 


- 2000 


103.5 _ 107.5 


1052 _ 1054.5 


X-ray emission 














zone 




15-50 


30 


- 100 


10^ - 10* 


1052.5 _ 1054.5 


(alternative) 


1012.5 


<1.2 




1000 


109.5 


1047.5 



Notes. — In the table, the model parameters R: source radius, T: bulk Lorentz factor, typical 
Lorentz factor of relativistic electrons, rig,: number density of emitting electrons, and A^e^ total 
number of electrons. 



Similarly, although order-of-magnitude variations in brightness are seen in the X-rays, 
the detail ed radio light curve does not reveal the co incident variations that would be expected 
for SSC (IZauderer et al.l l201ll : iLevan et al.l l201ll ). We therefore conclude that the X-ray 
emission must originate from a region separated than the radio emission region. 

While observations made in the first few days cannot discriminate between the two 
alternatives for the X-ray emitting region, we show in ^ below that the re-brightening 
observed at radio flux after tens of days can be naturally explained as resulting from material 
collected by the X-ray emitting plasma, provided that it travels at F 10 — 20. Thus, the 
first model, in which the X-ray emission region is separated from the radio emission region 
in velocity space is preferred. This leads us to a proposed two component jet, in which the 
fast. X-ray emitting material, is surrounded by a slower, radio emitting material. A cartoon 
demonstrating this model is shown in Figure [3l 



3. Late time radio evolution 



The two component jet model presented above has a distinct prediction. As the plasma 
expands into the interstellar material (ISM), it collects material, slows and cools. This sit- 
uation is similar to the afterglow phase in GRBs. Thus, using similar assumptions, one 
can predict the late time synchrotron emission from the decelerating plasma. One notable 
difference from GRB afterglow, though, is that the X-ray emitting plasma is mildly rela- 
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X-rays 




Fig. 3. — Sketch of geometrical configuration and emission regions for TDE Swift 1644+57. 
The inner jet has a Lorentz factor F ^ 15, and is responsible for the early times X-ray 
emission. The outer jet has Lorentz factor < 1.2, and is the source of the early times 
radio emission. As the inner jet propagates through the ISM it collects material and cools. 
Synchrotron emission from the collected material is the source of the re-brightening seen in 
radio frequencies after ^ 30 days. 



tivistic, while the radio emitting plasma is trans-relativist ic. Thus, when calculatiii g the 
dynamics, one cannot rely on the ultra- relativistic scheme (IBlandford fc McKedll976l ). but 
has to consider the transition to the Newtonian regime. 



3.1. Dynamics and radiation from an expanding jet 



The dynamics of plasma expanding through t h e ISM is well studied i r i the litera- 



ture (IBland 



Huang et al. 



brd fc McKe"3 Il976l: IXatz fc PiranI EqQtI: Ichiang fc Dermer 



1999[ Ivan Paradiis et allboool : 



Pe'ei 



20121 



1999; Piran 1999 



Nava et al.ll2012l ). Here we briefly 



review the basic theory, which we then utilize in calculating the expected late time radio 
emission from TDE Swift 1644+57. 



We assume that by the relevant times (tens of days), the reverse shock had crossed the 
plasma, and thus only the forward shock exists. The evolution of the bulk Lorentz factor of 
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the expanding plasma is given by (jPe'erl |2012| ) 



dT ^ 7(r^-i)-(7-i)r/3^ 

dm M + em + (1 - e) m [27! - (7 - 1) (1 + T-^)] ' ^ ^ 

Here, M is the mass of the ejected matter, m is the mass of the collected ISM, F is the bulk 
Lorentz factor of the flow, 7 is the adiabatic index and e is the fraction of the shock-generated 
thermal energy that is radiated (e = in the adiabatic case and e = 1 in the radiative case). 
Note that equation (I7j) holds for any value of F, both in the ultra-relativistic (F ^ 1) and 
the sub-relativistic (/3 ^ 1) limits. 

Under the assumption of constant ISM density, the collected ISM mass is related to the 
distance dR and observed time dt via 

dm = AuR'^nisM'mpdR; 

dR = F/3c (F + F/3) dt = j^cdt, ^ ^ 

where rup is the proton's rest mass and nisM is the ISM density. In deriving the second 
line in equation (jS]), we have explicitly assumed that the observed photons are emitted from 
a plasma that propagates towards the observer. A more comprehensive calculation which 
considers the integrated emission from diffe rent angles to the line of sight results in a similar 



solution, up to a numerical factor of a few ( IWaxmanI 1 1 9 9 Tt iPe'er &: Wijersll2006l ). 



In order to predict the synchrotron emission, one needs to calculate the magnetic 
field, B and characteristic electron's Lorentz factor, 7ci. This calculation is done as fol- 
l ows. By solving the shock jump conditions, one gets the energy density behind the shock 



(iBlandford fc McKed 119761 ) 



U2 = (T - 1) ^'^ ^^ riisMmpC^ . (9) 
7-1 

Equation (Q is exact for any velocity, including both the ultra-relativistic and the Newtonian 



limits. A useful approximation for the adiabatic index is 7 = (4F + 1)/(3F) (e.g., iDai et al. 
Il999l ll^. In the ultra-relativistic limit, F 3> 1, 7 = 4/3 and equation ([9]) takes the form 
U2 ~ 4F^nigMmpC^; While in the Newtonian limit, /3 ^ 1 and 7 = 5/3, equation becomes 

The shock-generated magnetic field assumes to carry a fraction ee of the post-shock 
thermal energy, B'^/Sn = eBU2, 



^ _ I (327reBnisM"^pc2)^/^F (relativistic limit), ^^^^ 
(167reB?^isM"^pC^)^^^ /3 (Newtonian limit). 



more accurate formula appears in lPe'eij (|2012l ). 
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Similarly, a constant fraction ee of the post-shock thermal energy is assumed to be carried 
by energetic electrons, resulting in 7eimeC^ = ee{u2/n2), where the number density in the 
shocked region is n2 = nisuil^ + l)/(7 ^ !)• This leads to 

{CeF ( ^ j (relativistic limit), 
ee^- f (Newtonian limit). 

In order to calculate the observed peak frequency and flux, we discriminate between two 
cases: 



V. 



p 

e.g., 



(i) In the optically thin emission, Va < i.e., t,^,,, < 1, the peak synchrotron frequency is 
= {^lATT){q B/mMV and the peak flux = F,^ ^ {Nj47rdl){2aTm,c^BT /9qe 



,ob ,,ob 



Sari et al.lll998l ). Here, di is the luminosity distance and A^e oc is the total number 
of radiating electrons (originating both at the explosion, as well as collected ISM). In the 
limit where the total number of swept-up ISM material is much larger than the original 
ejected material, analytic scaling laws can be obtained in the ultra-relativistic and Newto- 
nian limits. In the relativistic regime, F 3> 1, F oc t~^^^ and R oc t^^^, and thus 
and oc On the other hand, in the Newtonian limit, f3 oc t~^^^ and R oc t^^^, one finds 



oc t-^ and K oc t^/^ 



p 

(ii) In the optically thick regime, z/^ > z/^, i.e., r^^ > 1, the observed peak frequency Up 
and the peak flux F^^ are at the self absorption frequency, z/^. For a power law distribution 
of electrons above 7ei with power law index p, these are given by 

,,ob _ ,.ob _ ,.ob^2/(p+4) _ 3 9ei3^2p^2/(p+4) 

/ -E^ o -E^ (12) 

F— F (i^l 2 CTT"tcC RFt 

The optical depth at z/^, t^j^ is calculated as follows. We first point out that the comov- 
ing number d ensity of the electrons in the shocked plasma frame is n2 — 4FnisME Using equa- 



tion f6.5 2l inlRvbicki fc LightmanI ( 119791 ) . we find that r^,„ = a^^{R/T) = fip)nisMB-^%^^R 
(see also lPe'er fc Waxmanlbooi THere. f{p) = {p-l){87rV3qe/9)2P/^T{{2p + 2)/12)T{{3p + 
22)/12) is a function of the power law index, p of the accelerated electrons above 7ei, and 
F(x) is F function of argument x. Using this result in equation ( !T2|) . one can obtain analytic 
scaling laws for the temporal evolution of z/°^ and F^^ in both the ultra-relativistic and New- 
tonian limits, in a similar way to the derived expressions in the optically thin limit above. 
In the relativistic regime, z/°^ oc t-(3p+2)/2{p+4) ^ ^-2/3 ^ 2) and F^^ oc t5{i-p)/2{p+4) ^ 



^The equation relating the downstream (712) and upstream (rii) number densities, n2 = (7r + l)ni/(7— 1) 
holds for any strong shock, at any velocity, both in the relativistic and Newtonian regimes. 
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t-5/i2(^p = 2). In the Newtonian limit, one finds z/°^ oc t(2-3p)/{p+4) ^ f'^/^ {p = 2) and 

F,^ oc t(47-32p)/5(p+4) ^-17/30 (p = 2). 

3.2. Radiative cooling of the electrons 

In the calculation presented in §3.11 (in particular, equation ( ITT]) above), we considered 
heating of the ISM plasma as it crosses the shock front. Once the ISM particles cross the 
shock front, they lose their energy by radiative cooling. Hence, their typical Lorentz factor 
7ei decreases with time. Similarly, the original ejected material from the tidally disrupted 
star cools with time. 

Within the context of our double jet model, this cooling has a more pronounced effect 
on the lightcurve originating from the outer jet, that is responsible for the early time radio 
emission. This is because this region is composed of a dense plasma, that propagates at trans- 
relativistic velocities (see Table [T]). Thus, the contribution of the swept-up ISM material 
to the emission during the first ~ month from this region is minor, as opposed to the 
contribution of the collected ISM to the emission from the inner jet. This can be seen by 
noting that during Newtonian expansion, R oc t"^^^. Thus, between five and thirty days, 
the radius of the radio emitting region is increased by a factor ^ 2, and therefore cannot 
exceed ^ 2 x 10^^ cm. As a result, unless the ISM density is much larger than ~ lO^cm"^, 
contribution from the collected ISM to the emission from the outer jet is subdominant. 

The radio emitting electrons cool due to synchrotron and inverse Compton scattering. 
The radiated power from a relativistic electron with Lorentz factor 7ei is P = Pgyn + -Pic = 
(4/3)ccrr {Ub + Uph) 7ei, where Ub and t/ph are the magnetic and radiative field energy den- 
sities, respectively. 

Both Ub and f/ph decrease with radius (and time). As discussed in §3.11 above, during 
the Newtonian expansion phase -B oc /3 oc t~^^^ (see equation f lTU]) ). and thus Ub oc t'^^^. 
Similarly, f/ph = L/{AttR^T^c), where the luminosity is L oc N^-f^B'^T^ oc oc t"^^/^ 
(assuming that most radiative particles are from the original ejecta). This leads to f/ph oc 
^-22/5^ namely the IC component decreases much faster than the synchrotron component. 

We can therefore write P ~ (4/3)c(TTf/B,o(^Ao)^^^^7ei) where Uba is the magnetic field 
energy density at fiducial time, to which is taken in the calculations below to be five days 
(observed time). The characteristic particles Lorentz factor at any given time t > to is given 
by 
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where 7ei,o is the electron's Lorentz factor at to- The temporal evolution of the electron's 
Lorentz factor and the decay law of the magnetic field allow calculation of the late time 
evolution of both the peak radio frequency and peak radio flux from the outer jet region. 



4. Interpretation of the late time radio emission of Sw J 1644+57 

The temporal evolution of the radio flux and peak radio frequency of Sw J1 644+57 up to 



216 d ays from the initial outburst are presented in Figures|ll[5](data taken from lBerger et al. 



20121 ). Three separate regimes are identified in both figures: (I) At early times, < 30 days, 
the flux decreases from ~ 40 mJy at ~ 5 days to ~ 15 mJy at ~ 30 days. During this period, 
vf" rapidly decays, with a decay law consistent with oc and a ~ 2.0. (II) Between 
30 — 100 days, the flux increases by a factor of ~ 1.8 (from ~ 20 mJy to ~ 35 mJy). At the 
beginning of this epoch, at ~ 30 — 45 days, the radio peak frequency p'^ increases by a 
factor of a few, while during the rest of this epoch it shows a similar decay law as is seen at 
early times, v'^ oc with a ~ 2.0. (Ill) Finally, at very late times, > 100 days, the flux 
decays again. This decay is accompanied by a slow decay in the peak radio frequency, much 
slower than the decay observed at earlier phases. W e point out that the increase in the radio 



flux observed at epoch (II) led iBerger et al.l (120121 ) to suggest that late time central engine 
activity may take place. 

Based on the discussion presented in the previous sections, we suggest the following 
interpretation to the late time radio emission of Sw J1644+57: (I) During the early decay 
phase, < 30 days, the emission is dominated by synchrotron radiation from the same 
electrons that emitted the radio emission at early times; this is the outer jet component in 
our model (see Figure [3]). This emission is thus a continuation of the emission observed at 
early times. The observed decay of both the flux and peak frequency is due to the decreasing 
of magnetic field and the radiative cooling of these electrons. (II) At ~ 30 days, there is a 
transition: the inner jet plasma, that originally emitted the X-ray photons, expands into the 
ISM, collects ISM material and cools. At this time, synchrotron emission from this plasma 
becomes the dominant component at radio frequencies. Initially, the inner jet propagates 
at relativistic speeds, with F > 15 (see discussion in §2.21 and Table [1]). However, as it 
propagates into the ISM, the plasma collects material from the surrounding ISM and slows; 
the collected material contributes to the radio emission, resulting in an increase in the radio 
flux. (Ill) Finally, at ~ 100 days, the emission becomes optically thick (self absorption 
frequency is larger than the peak frequency), which causes the late time decay. 

Fits to the temporal evolution of the peak radio frequency and peak flux are shown in 
Figures H] and O In producing these fits, we use the following parameters, which match the 



- 15 - 



T 1 — I — I — I — r-| 1 1 1 1 1 — I — I — i—p 




'I I I I I I I I I I I I I I I I i_ 

10' 10^ 



Observed time [days] 



Fig. 4. — Temporal evolution o f peak flux of the radio emission of Swift 1644+57. Data is 
taken from iBerger et al.l ( 120121 ) . Dash-dotted line is the contribution of the outer (slower) 
jet, while solid and dotted lines are the contribution from the inner (faster) jet. As the 
inner jet propagates into the ISM, it collects material and cools; thus, after ~ 30 days, radio 
emission from this region dominates the radio flux. At ~ 100 days, the peak flux enters the 
optically thick regime (r^.^ > 1), which causes the decay see at these times. The values of 
the free model parameters used are presented in the text. 



early and late time properties of the flow. For the outer jet, we use initial expansion radius 
(at observed time = 5 days) Rq = 10^^ cm. The typical electron's Lorentz factor is taken 
to be 7ei_o = 190, magnetic fleld Bq = 4.5 G, and initial bulk Lorentz factor r(tQ'^) ~ 1.18. 
These values correspond to ec ~ 0.58 and cb ~ 0.13, and imply total number of radiating 
particles Ne{t()) = 1.25 x 10^^. These values are consistent with the flndings in §2.11 (see 
Tabled]). These parameters result in initial peak synchrotron frequency z/°q = 8.0 x 10^^ Hz, 
and peak observed flux F^^^q = 39.4 mJy. 

For the inner jet, we flnd that the best flt is obtained when using initial bulk Lorentz 
factor r(tg^) = 17, = 0.54, eg = 0.54 and r(tg^) ~ 1.2 x 10^^ cm. We point out that the 
value of is not constrained by early X-ray data. Moreover, here the collected ISM plays 
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Fig. 5. — Temporal evolut ion of peak f r equen cy for the radio emission of Swift 1644+57. 
Observed data taken from iBerger et al.l (120121 ). Lines have the same meaning as in Figure 
m Three distinctive regimes are clearly seen. 



a significant role. Thus, we could constrain the ISM density to be riisu = 2.7 cm~'^. We 
further assume that the shocked ISM have a power law distribution above 7ei with power 
law index p = 2. When calculating the dynamics in equation (j7]), we consider adiabatic 
expansion, namely e = 0. 

As the inner jet propagates into the ISM it collects material and slows. Using Equations 
[71 El we find that at observed time = 100 days, it reaches radius R ~ 6.0 x 10^'' cm, 
with bulk Lorentz factor F ~ 1.15. At this time, the magnetic field is B ~ 0.19 G and 
the characteristic electrons Lorentz factor is 7ei ~ 150. This result in z/°^ = 2.06 x 10^" Hz 
and F^^ = 34.4 mJy. However, at this time, t^,^ ~ 1, which implies that from this time 
onward synchrotron radiation is in the optically thick regime, and its late time properties 
are described by equations (fT2!) . While the dotted line in in Figures HI [5] presents the model 
fits in the optically thin regime, the solid line in these figures represent the expected peak 
frequency and flux in the optically thin region at early times, and optically thick regions at 
> 100 days. 
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5. Conclusions and Discussions 

In this paper, we studied the early and late time emission of the TDE Sw J1644+57, 
both at radio and X-ray frequencies. Based on our analysis, we propose a two-component jet 
model (see Figure |3]) that fits the observations. Our model explains, for the first time, both 
the early time X-ray and radio emission, as well as the complex late time radio behavior. 

By analyzing the early time {t°^ < 5 days) data, we conclude that the radio and X-ray 
photons must have separate origins (see ^ Table [1]). We are able to put strong constraints 
on the properties of the radio emitting plasma (see §2. H and Figured]), and somewhat weaker 
constraints on the properties of the X-ray emitting plasma ( §2.2l and Figure^]). Stronger con- 
straints on the initial properties of the X-ray emitting plasma are obtained when considering 
late time radio emission. 

Our jet within a jet model naturally explains the complex temporal evolution of the 
radio emission at late times (up to > 200 days). We solved in [12] the dynamics and expected 
radiation from the jet propagating into the ISM. We stress that as the initial Lorentz factor 
of the inner jet is mild (our best fit gives F ^ 17), one has to consider the transition between 
the relativistic and the Newtonian expansion phases, and cannot rely on analysis of ultra- 
relativistic outflows, as is the case in GRBs. We further consider the radiative cooling of 
particles behind the shock front. 

We demonstrated how our model can be used to fit the data in §1] Our key idea is that 
the observed signal can be split into three separate regions: at early times, radio emission 
is dominated by the outer jet. The decay of the peak frequency and flux is attributed to 
radiative cooling of the electrons and the declining of magnetic field. Between 30 — 100 days, 
the radio emission is dominated by the inner jet, in the optically thin regime. The inner 
jet propagates at relativistic speeds and collects material from the surrounding ISM. The 
addition of the collected material results in the increase of radio flux. Finally, at very late 
time, t°^ > 100 days, the radio emission is dominated by the inner jet, in the optically thick 
regime. This causes the observed late time decay of peak flux and frequency. 

When analyzing the early time data, we use the common interpretation of the radio 
emission as originating from synchrotron radiation. We find that the bulk Lorentz factor of 
the synchrotron emitt ing plasma is F < 1.2. Our analysis method is similar to that used by 



Zauderer et al.l ( 120111 ). albeit being more general, as we avoid the equipartition assumption 
used in that work. The X-ray emission is interpreted as IC scattering of the radio photons, 
however our analysis indicates that the bulk motion of the X-ray emitting region is much 
higher, F > 15; alternatively, the emission radius is much smaller, r ~ 10^^'^ cm. We thus 
conclude that radio and X-ray emission zones are completely separated, either spatially or 
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in velocity space. This is consistent with the separate temporal behavior observed at early 
times, as the observed X-rays maintained a more constant level after the first 48 hours (albeit 
with episodic brig htening; and fading spanning more than an order of magnitude in flux; see 
Levan et al.ll201l[ ). whereas the low frequency emission decreased markedly. 



Our calculated X-ray luminosity can match the observed X-ray luminosity of TDE Sw 
J1644+57, Lx ~ 10^^'^ erg sec~^, exceeding the radio emission by a factor ~ 10^ — 10''. If the 
emission is isotropic, the X-ray luminosity of Sw J1644-I-57 corresponds to the Eddington 
luminosity of an ~ 10^ Mq MBH, whi ch is incompatible with the upper limit ~ 10^ Mq of the 
MBH mass derived from variability (IBloom et al.l 1201 ll : iBurrows et al.l 120111 ) , so the source 
is required to be relativistically beamed. We have obtained the bulk Lorentz factor of the 
relativistic jetted outburst F ~ 15, which is very close to the typical value inferred in blazars. 
Therefore, if the beaming angle of the jet 6j l/F ~ 0.1, the beaming-corrected luminosity 



X 



10 



45 



erg sec 
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cosOj) is the beami ng factor) becomes consistent with the 



Eddington luminosity of a ~ 10^ Mq MBH (see also iBloom et al.l 1201 ll ). 



Berger et al.l ( 12012| ) present the continued radio observations of the TD E Sw J1644-I-57 



extend ing to ~ 216 days. In their work, they fitted the data using the model of lGranot fc Sari 
( 120021 ). and concluded that the re-brightening seen after ~ 30 days cannot be explained in 
the framework of that model. They thus conclude that an increase in the energy by about an 
order of magnitude is required. As discussed above, the re-brightening is very natural in our 
scenario due to the increase in the collected material when the jetted outburst propagates 
through the ISM. 

The idea of jet within a jet was suggested in the past as a way to explain the morphology 
of high energy emission from active galactic nuclei (AGNs) ( iGhisellini et al.ll2005t iHardcastle 
20061 : iJester et al.ll2006t ISiemiginowska et al.ll2007n. and in the context of GRBs in explaining 



the break observe d in the afterglow light curve (IRacusin et al.ll2008l : Ide Pasauale et al. 



Filgas et al.ll201l[) . as well as some of the properties of the prompt emission (iLundman et al. 



20091: 



2OI2I ). While the theory of jet launching is still incomplete, clearly, jets, being collimated 
outflows are expected have lateral velocity gradient. The analysis done here suggests that 
such a velocity gradient exists in jets originating from TDEs. In principle, the outer jet may 
also represent a cylindrical boundary layer owing to the interaction of the inner jet with the 
ambient gas. 

While Sw J1644+57 is the first TDE event from which the existence of relativistic jets 
is inferred, and the most widely discu ssed one in this co ntext, it is possible that other 
such events were detected. Recently, ICenko et al.l (120121 ) reported a second event, Sw 
J2058. 4+0516 which is a potential candidate for a relativistic fiare. While currently, no 
late time radio lightcurve is currently available, we can predict that if such a lightcurve 
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becomes available it should show the same complex behavior of Sw J1644+57. 
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